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ABSTRACT 



This paper introduces an important new canonical set of functions for solving 
Lanchester~type equations of modern warfare for combat between two homogeneous forces 
with power attrition-rate coefficients with **no offset.*’ Tabulations of these func- 
tions, which we call Lanchester-Clif ford-Schlaf li (or LCS) functions, allow one to 
study this particular variable-coefficient model almost as easily and thoroughly as 
Lanchester’s classic constant-coefficient one. The availability of such tables is 
pointed out. We show that our choice of LCS functions allows one ‘to obtain important 
information (in particular, force-annihilation prediction) without having to spend the 
time and effort to compute force-level trajectories. Furthermore,^ we show from theoret- 
ical considerations that our choice is the best for this purpose. These new theoreti- 
cal considerations apply in general to Lanchester-type equations of modern warfare and 
provide guidance for developing other canonical Manchester functions (i.e. canonical 
functions for other attrition-rate coefficients). Moreover, our new LCS functions 
provide valuable information about various related variable-coefficient models. Also, 
we introduce an important transformation of the battle’s time scale that not only many 
times simplifies the force-level equations but also shows that relative fire effective- 
ness and intensity of combat are the only two weapon-system parameters determining 
the course of such variable-coefficient Lanchester-type combat. 



0. Introduction 



In an earlier paper ^ we developed some elements of a mathematical theory for 
solving variable-coefficient Lanchester-type equations of modern warfare for combat 
between two homogeneous forces and introduced canonical hyperbolic-like Lanchester 
functions for constructing their solution. Unfortunately, with only these previous 
results one is limited to computing (however conveniently it may be done)- force-level 
trajectories and cannot gain a real understanding of qualitative model behavior (e.g. 
force annihilation) without the excessive labor of extensive numerical computations (and 
only then for specific values of model parameters). Since our earlier work important 
mathematical discoveries have been made about* the qualitative behavior of this combat 
model, and we wish to show how these new results allow one to parametrically analyze 
combat modelled by power attrition-rate coefficients ( see Section 1 below) with some- 
what the same facility as he cai\ study F. W. LANCHESTER ’ s classic constant-coeffi- 
cient model. In order to obtain this analysis capability, however, one must redefine 

the Lanchester-Clif ford-Schlaf li (or LCS) functions, which we introduced in reference 24. 

[241 

In our earlier paper ( see TAYLOR and BROWN ) we gave various examples of 

hyperbolic-like Lanchester functions (in particular, the LCS functions, which arise 

from power attrition-rate coefficients with ”no offset’*). Subsequent research by 

F271 

TAYLOR and COMSTOCK has revealed, however, that these canonical LCS functions must 
be redefined to permit force-annihilation prediction from initial conditions without 
having to spend the time and effort to compute force-level trajectories. It 
then became obvious that the entire topic. of representii^ the solution to such Lan- 
chester-type equations in terms of general Lanchester functions (GLF) should be 
critically reexamined. Thus, the purpose of this paper is to present new general con- 
siderations for the selection of canonical Lanchester functions and then apply this 
theory to the special case of power attrition-rate coefficients with **no offset” 
(modelling, for example, weapon systems with the same maximum effective range) to obtain 
new LCS functions. With the availability of tabulations of these new LCS functions 



( see Section 6 below) , one can study this model almost as easily and thoroughly as 
Lanchester^s classic constant-coefficient one. 

These power Lanchester (i.e. LCS) functions are significant not only because 
they correspond to attrition-rate coefficients modelling a large class of combat 
situations of interest but also because they yield valuable information about other 
related canonical Lanchester functions, e.g. the offset power Lanchester functions 
(see Note 1 and Section 8 below). This information is, of course, equivalent to know- 
ledge about model behavior (e.g. force annihilation). As a result of our work here 
one can parametrically analyze variable-coefficient Lanchester-type models for combat 
between two homogeneous forces with somewhat the same facility as the classic con- 
stant-coefficient one. Such models are important for developing insights into the 
dynamics of combat (see BONDER and TAYLOR ^ ^ , and Section 1 below). 

The organization of this paper is as follows. First, we present the variable- 

coefficient Lanchester-type model that we study in this paper. Next, we discuss the 

representation of the time history of the force levels for this model in terms of 

general Lanchester functions (GLF) . We show that there are essentially only two kinds 

of GLF, (I) exponential-like GLF and (IT) hyperbolic-like GLF, and that the former (I) 

provide essential force-annihilation-prediction information about the latter (II). 

Then we explain why we have chosen to use the hyperbolic-like GLF to construct the 

model’s solution and why the power Lanchester (or LCS) functions introduced by Taylor 
[241 

and Brown must be redefined. Next, we show how the analysis of, for example, the 
X force-level equation is simplified by transforming the independent variable t to 
normalize the battle’s time scale by the intensity of combat. We then introduce our 
new definition of Lanchester-Clif f ord-Schlaf li functions and show how they arise in 
solving the transformed X force-level equation. Availability of tabulations of 
these new LCS functions is discussed, and some uses of the tabulations are illustrated. 
Finally, insights gained into the dynamics of combat between two homogeneous forces 
from these developments are discussed. 
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1. Variable-Coefficient Lanchester-Type Equations of Modern Warfare 

In this paper we consider the following idealized model for combat between two 
homogeneous forces ( see Note 2) 

dx/dt = -a(t)y, dy/dt = -b(t)x, (1) 

with initial conditions 

x(t=0) = Xq, and y(t=0) = y^, 

where t=0 denotes the time at which the battle begins, x(t) and y(t) denote the 

numbers of X and Y at time t, and a(tX and b(t) denote time-dependent Lan- 

• • 

Chester attrition-rate coefficients. We will refer to (1) as variable-coefficient 

Lanchester-type equations of modern warfare in honor of the pioneering military opera- 

roil f9Ai 

tions research work of F. W. Lanchester^ ^ ( see TAYLOR^ and Taylor and Brown^ 

Other forms of Lanchester-type equations appear in the literature, but we will not con- 
sider these here ( see DOLANSKY^^^^ and Taylor The Lanchester-type equations (1) 
yield the X force-level equation 

d^x/dt^ - {d £n a(t)/dt}dx/dt - a(t)b(t)x = 0, (2) 

with initial conditions 

^ x(t=0) = Xq, and { [l/a(t) ]dx/dt = -y^. 

Although combat between two military forces is a complex random process, such 

an idealized deterministic model of the combat attrition process is frequently employed 

[9] 

to provide insights into the dynamics of combat ( see , for example, BONDER and FARRELL 

Bonder and Honig^^^^, TAYLOR and PARRY^^®^, or WEISS^^^.^). The reader may consider (1) 

to model combat in which both sides use aimed fire and target acquisition times are inde 

pendent of the numbers of firers and targets ( see Note 3) . New operations research 

roi fi21 

techniques ( see , for example. Bonder and Farrell^ \ and CLARK ) for forecastxng 
temporal variations in fire effectiveness (caused by, for example, changes in force 
separation, combatant postures, target acquisition rates, etc.) have generated interest 
in such variable-coefficient combat formulations. 
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Without loss of generality, we may take a(t) and b(t) to be of the form 



a(t) = k^g(t), and b(t) = k^h(t), (3) 

where g(t) and h(t) denote the time-varying factors of a(t) and b(t) such that 
a(t)/b(t) = g(t) = h(t). In other words, k^ and k^ denote "scale” 

factors chosen so that the case of constant coefficients corresponds to g(t) = h(t) = 

A large class of tactical situations of interest can be modelled with the follow 
ing general power attrition-rate coefficients 

a(t) = k^(t+Kg)^, and. b(t) = k^Ct+Kg+K^)'', (4) 

where Kg,K^ ^ 0. The modelling roles of and are discussed in Taylor and 

[241 

Brown . We will call K the starting parameter , since it allows us to model 
(with p,v ^ 0) battles which begin within the maximum effective ranges of the two 
systems. We will call the offset parameter , since it allows us to model (again, 

with y,v ^ 0) battles between weapon systems with different effective ranges. Restric 
tions that must be placed in y and v, which are not necessarily integers, are 
discussed below. 

Let us take a few moments to motivate our above notation and further indicate 
possible applications of our results. Consider BONDER’s^^*^^ constant-speed attack 
on a static definsive position modelled by 

dx/dt = -a(r)y = ~aQ(l-r/R^)^y , dy/dt = -B(r)x - -8Q(l-r/Rg)^x, (5) 

where y,v ^ 0 and R denotes the maximum effective range of the Y weapon system. 

a 

Then the starting parameter and the offset parameter are given by i 

Kg = (R^-Rq)/v, and 

where R^ denotes the battle's opening range and v > 0 denotes the constant attack 
speed. Hence, K^,K« ^ 0 o ^ R ^ R^. By considering (6) and Figure 1, the reader 

w D p Ot U 

should have no trouble in understanding our terminology for Kg and K^. In the model 
(5) p, for example, is used to model the range dependence of Y's attrition-rate 
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Figure !• Explanation of offset parameter and starting parameter for power attrition-rate 

coefficients modelling constant-speed attack. [Notes: 1. The maximum effective ranges of 

the two weapon systems are denoted as R and R. . 2. The opening range of battle is 

ot p 

denoted as R- and (as shown) < minimum (R ,R.). 3. The offset parameter is given by 

U U CX p 

Kp, =■ (Rq-R )/v. 4. The starting parameter is given by K„ = (R -R»)/v.] 



coefficient ( see Figure 2). Observing that range is related to time by r(t) = R^-vt, 

we readily see that the longest the battle can last is given by t = at which 

max U 

time zero force separation is reached. 

When the offset parameter is equal to zero (i.e. Rq “ 0) > then the coefficients 
(4) reduce to 

a(t) = k^(t+Kg)'^, and b(t) = kj^(t+Kg)'’. (7) 

We will refer to (7) as power attrition-rate coefficients with **no offset .** The purpose 

[241 

of this paper is to extend our previous results and introduce new power Lanchester 
functions that allow more information to be more conveniently extracted from the model 
(1) with coefficients (7). Specifically, one would want to obtain information such as: 
(Ql) Who will ’’win**? Be annihilated? 

(Q2) How do force levels decrease over time and how many survivors will the 
winner have? 

(Q3) How do changes in the initial force levels and weapon system parameters 
affect the outcome? Is concentration of forces a good tactic? 

(Q4) How long will the battle last? 

To conveniently answer questions (Ql) , (Q3) , and (Q4) a redefinition of the Lanchester- 
Clif ford-Schlaf li (or LCS) functions is required. Moreover, not only are results for 
the coefficients (7) of interest in their own right but they also provide much valuable 
information about the general case (4). 

2. Representation of Solution in Terms of General Lanchester Functions 

In this section we discuss how to construct the solution to the X force-level 

equation (2) in terms of certain basic building blocks that we have chosen to call 

general Lanchester functions (GLF) . We feel that these GLF should be chosen according 

to the guidelines shown in Table I. Special cases of these general considerations have 

been given by Taylor and Brown^^^^ and Taylor and Comstock^^^ ^ . 

T T 

Let us introduce the GLF x = (x- x^) and y = (y^ y^) which satisfy 

^ 1 Z ^ 1 z 

x= ka(t)L;^, £= (l/k)b(t)Lx, (8) 
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(range). 3. The opening range of battle is denoted as = 1250 meters 



TABLE I. Requirements for General Lanchester Functions 



(Rl) : They can be used to construct the solutions to the X and Y 

force-level equations • 

(R2): They should be as ’’simple” as possible. 

(R3) : A given set of functions should apply to as large a class of 

battles as possible. 

(R4) : They should be nonnegative. 

(R5) : They should reduce to elementary transcendental functions in 

special cases such as a constant ratio of attrition-rate 
coefficients . 

(R6) : They should provide as much information as possible about model 

behavior. 
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where without further specification k may be any positive constant and L is any 
2x2 square matrix such that = !• [The initial conditions for (8) are to be 
chosen so that the requirements of Table I are met and are discussed below.] Then x 
satisfies the vector equation 

x-{d£na(t)/dt}x-a(t)b(t)x = 0. (9) 

In other words, both x^ and satisfy the X- force-'level equation (2), and 

similarly y^^ and y^ satisfy the Y force-level equation. 

Let us now investigate all the possible forms for the above GLF. Intuitively, 
we would expect two possibilities (keeping the requirements of Table I in mind) : 
exponential-like functions (one strictly increasing and the other strictly decreasing) 
and hyperbolic-like functions. [We are reminded of these two possibilities by the 
well-known constant-coefficient results.] Two such types of GLF appear in Taylor and 
Brown ^ (hyperbolic-like functions) and in Taylor and Comstock^ (exponential-like 
functions) . We will show that these are the only two possibilities if the requirements 
of Table I are to met and show the relationship between these two types of GLF. 

It is readily shown that any 2x2 matrix such that L^ = I must take one of 
five forms. 



LEMMA 1: If = I, 

/a S \ 


■ then L 


must 

/ 


take one 
1 0\ 


of the 

i 


following 
r-l Ov 


five 


forms : 


(A) 


I ) with 

V[l-a2]/B -a/ 


B 0, 


(B) ( 


, J- 


(C) j 


.Y J’ 


(D) 


L = I, 


or 



(E) L = -I, where a, S, and y are unrestricted with the exception that 



& ¥ 0 . 



T 

If in addition L = L and |l| = -1, then L is an orthogonal matrix and must be 
of the form 

( cos (() sin (J)\ 

). ( 10 ) 

sin (J) -cos <{)/ 

It seems reasonable to give our Lanchester functions symmetry by requiring that L = 

[241 

We observe that the hyperbolic-like GLF of Taylor and Brown correspond to 
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In order that (R4) of Table I be met for a constant ratio of attrition- 






rate coefficients, i.e. 



a(t) = k h(t), 

a 



and b(t) = kj^h(t) , 

r o / 1 

we know ( see Taylor and Brown^ that we should choose k = /k^/k^ so that 

k = a(t)L;^, i = Ajk^ b(t)Lx. 

The general solution to the X force-level equation (2) is given by 



x(t) = C X. 



( 11 ) 



( 12 ) 



(13) 



Introducing the 2^2 matrix 



S(t) 



T 

X (t) 






(lA) 



^;^‘(t)L 

we obtain from the initial conditions to (2) that 



- 1 , 



C = S (t=0) 



(A ) 



(15) 



Observing that 



^ |S(t)l = 



x^t) 

X^t)L^^ 



/(t) 

•T/ 

X (t)L 



we readily see that |s(t)| = constant Vte(tQ,+“), where |s| 
of the square matrix S. Thus, we may take 

|s(t)| = constant . 



denotes the determinant 



(16) 



Let us also observe that 



A^/k^ a(t) |s(t) I = W(Xj^,X2) 



(17) 



where W(Xj^,X 2 ) ^®"°tes the Wronskian of x^^ and X 2 < 



Now let us subject the fundamental system of solutions x to linear trans- 
formation 

X = Ax, (18) 

such that the form of the equations (12) remains invariant, i.e. 
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(19) 



£= A^/k^ b(t)Lx, 

^here L again is such that L =1. If L is given, it follows that 

£ = LAL;^. (20) 

furthermore 

S^(t) = AS^(t), (21) 

50 that 

|s(t)l = |A|*|s(t)| = constant. (22) 

le also observe that W(x^,X 2 ) = [a|w(x^,X 2 ) Considering the quotient of the two 
general Lanchester X-functions (GLXF) 

n(t) = x^/x^, (23) 

7e see that under the linear transformation (18) we have 

dn/dt = { I a| 

We now show that the only possible GLF that satisfy the requirements of Table I 
r T 

.corresponding to L = L ) are the exponential ones shown in Table II and the hyper- 
bolic ones shown in Table III. According to Lemma 1 if = I, then L must take 

me of five forms. It is impossible to have L = -I and satisfy (R4) of Table I ( see 
lote 4). If L = I and we try to specify the ’’simplest” initial conditions [i.e. 
ipecify initial conditions such that the GLF take the ’’simplest” form (satisfy (R2) of 
'able I)], we find that we may take L to have one of the three remaining forms (see 
fote 5). If we require that L be symmetric for simplicity [requirement (R2)], then 
4 is an orthogonal matrix with the form (10) . If in (10) we take cos (j) and sin <{> 
iqual to -1, 0, or 1 in order that the GLF take the ’’simplest” form, we find that 
:he only two distinct possibilities for L are L^ and L^ as given in Tables II 
md III ( see argument given in Note 5). Thus, we have shown that if we wish to con- 
struct the solution (13) to (2) by using GLF with the properties given in Table I, 
:here are essentially only two possibilities: the exponential-like GLF introduced by 
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Table II. Exponential-Like General Lanchester Functions 



" (E^(t;Q*) E”(t;Q*-)), " (EY(t;Q*) E^(t;Q*)) 



2. ^(t=tQ> =.(1/Q* 1), X^(t=tQ> = (1 Q*) 




4. - -2 
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Table III. Hyperbolic-Like General Lanchester Functions 



1. ^(t) = (S^(t) C^(t)), *;s^(t) = (S^(t) 



2. ^<t=to) = (0 1). = (0 1) 




4. |Sj^(t)|=-l 



c^Ct)) 
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Taylor and Comstock^^^^ and the hyperbolic-like GLF introduced by Taylor and Brown^^^^, 
We call the quantity Q*, which appears in the initial conditions for the 
exponential-like GLF, the parity condition parameter . It is chosen so that ( see Note 6) 

Ex(t;Q*), E^(t;Q*) > 0 for all finite t ^ t^. (25) 

It may be considered to be the enemy force equivalent of a friendly X force of unit 

[27] 

strength. Taylor and Comstock show how knowledge of the parity condition parameter 

allows one to predict force annihilation from initial conditions without explicitly 

computing force-level trajectories. We observe that the exponential-like GLF cannot 

f 27] 

be computed until one has solved the associated auxiliary parity-condition problem 

(i.e. knows how to predict force-annihilation). For this reason and others ( see Taylor 
[27] 

and Comstock ), the exponential-like GLF are mainly of theoretical importance. 
Moreover, in the next paragraph we show how the exponential -like GLF provide valuable 
force-annihilation information about the hyperbolic-like GLF. 

We now show that the limiting value of the quotient of the two hyperbolic-like 
GLXF, equal to the reciprocal of the parity condition parameter, 

i.e. (30) holds. We know that the two types of GLF are related by a linear transformatioi 



2h " 



(26) 



From (21) at t = 0, we have 



(“ n - A ^ ) . 

'10' W -Q*/ 



whence 

/1/2 -1/(2Q*)\ 

A = 

'Q*/2 1/2 ' 

Considering (23) and (26) , we see that 



(27) 



rijj(t) = {aj^j^rij.(t;Q*)+aj^2^/t^21^E^‘^’^*^'*’^22^* 



( 28 ) 



Recalling that ilc.(t;Q*) = E„(t;Q*)/E3(t;Q*) and lim E„(t;Q*) = 0, we see that 
L XX , X 
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( 29 ) 



so that by (27) and (28) 



lim n„(t;Q*) = 
t-H« ^ 



lim hj^(t) = 1/Q*. (30) 

t->+® 

rhis result (30) is highly significant because it leads to a computation method for 
determining Q*. Moreover, in the future we will show how the LCS functions introduced 
in this paper ( see Section 5 below) play a crucial role in such numerical determinations. 
Let us further note that (24) reads 

drijj/dt = {2/(Q*ng+l)^}dng/dt, (31) 

30 that both n„(t) and n„(t) are strictly increasing functions of t, since 

£j n 

3ng/dt = 2/k^/k^ a(t)/{E~(t;Q*)}2. 

3. Considerations for the Choice of General Lanchester Functions 

In Table I we give the general requirements that we feel should be placed upon 

jLF. These requirements are motivated by the properties possessed by the functions 

(namely, the exponential and hyperbolic ones) that one uses to construct the solution 

to (1) in the constant-coefficient case. 

We specify (R3) so that as few tabulations of GLF as possible will be required. 

X Y 

[Consequently, we specify the initial conditions for the GLF at tg = max(tQ,tQ), 

X 

^here t^ denotes the largest finite singular point on the t-axis for the X force- 
Level equation (2) ( see reference 24). Thus, a(t) and b(t) are positive continuous 
cunctions Vtc(tQ,+®). Since at most one of x(t) and y(t) can vanish in [0,-H») 

[ see Note 7), we have what the mathematician calls a nonoscilliatory solution to (1). 

[n this case we can construct the solution to (2) out of nonnegative components and 
/ill find it convenient to do so. 

As we have shown above, there are essentially only two types of GLF that 
satisfy the requirements of Table I: the exponential-like GLF of Table II and the 

lyperbolic-like GLF of Table III. We feel, however, that the hyperbolic-like functions 
ire to be preferred for two reasons: (1) they apparently are more convenient for 
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parametric studies in which one might, for example, want to vary initial force levels 

[241 

or some measure of relative fire effectiveness ; and (2) accurate values of the 

exponential-like GLF are, in general, difficult (in fact, essentially impossible for 

large values of t) to determine, since their initial conditions depend on the parity 

condition parameter Q* ( see Note 8). In terms of the hyperbolic-like GLF, the 

[241 

solution to (2) is given by 

x(t) = XQ{CY(t=0)C^(t)-SY(t=0)S^(t)}-yQ^iT7^ {Cjj(t=0)S^(t)-S^(t=0)Cj^(t) }. (32) 

We observe that for t^ < 0, for example, C^(t=0) > 1 and S^(t=0) >0 so that 

except for the quasi-autonomous case in which a(t)/b(t) = constant ( see Note 9), the 

[24] 

solution (32) only simplifies when t^ = 0 ( see Theorem 1 of Taylor and Brown ) 

( see Note 10) • 

Unfortunately, the power Lanchester (or LCS) functions introduced by Taylor and 
[241 

Brown were inappropriately defined to yield all the information sought about the 

combat model (1) with power attrition-rate coefficients (7) • In particular, the time 

• 

at which a side will be annihilated cannot be determined (without the explicit calcu- 
lation of the entire force-level trajectories) from the initial conditions, Subse- 

[2 

quent work by Taylor and Comstock has yielded a theory of force-annihilation prediction 
The purpose of the paper at hand is to redefine the power Lanchester functions in light 
of these subsequent results. We also thought it important to present the general 
considerations behind this selection of canonical Lanchester functions. 

Moreover, the form of the LCS functions is simplified and insight gained into 
the dynamics of combat by transforming the battle’s time scale. Thus, certain trans- 
formations of variables may be desirable in the development of hyperbolic-like GLF, 
and the specifications of Table III should be interpreted as being ’’symbolic” and not 
taken literally. 



16 



A Transformation to Normalize the Battle's Time Scale 
BY the Intensity of Combat 

ft 

Let . . . ds denote an indefinite integral, denote the relative effectiveness 
as R(t) , i,e, 

R(t) = b(t)/a(t), (33) 

and let K be an arbitrary constant to be conveniently chosen. Then Theorem 2 of 
[24] 

Taylor and Brown may be stated as follows. 

THEOREM 1: A necessary and sufficient condition to be able to transform 

the X force-level equation (2) by a transformation of the independent 
variable t into a linear second order ordinary differential equation with 
constant coefficients is that 

[d{An R(t)}/dt]//a(t)b(t) = CONSTANT, 



and then the desired substitution is given by 

rt 



T = K 



/a(s)b (s) ds. 



We observe that Theorem 1 says that we can transform the X force-level equation to 
a constant coefficient one if and only if d{^^n R(t)}/dx = CONSTANT. We also assume 
that the following condition holds. 



CONDITION (A) : 



a(s)ds and b(s)ds 



T 



t i t„. 



are bounded for all finite 



If Condition (A) is to hold, then for the power attrition-rate coefficients (7) we 
must have p,v > -1. 

Motivated by both Theorem 1 and the well-known constant-coefficient results, 
we introduce the new independent variable t defined by 



T 



/a(s)b(s) ds. 



(34) 
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By Condition (A) and the Cauchy-Schwarz inequality for integrals ( see p. 123 of BELLMAN 
the integral in (34) is well defined (i,e. bounded)* The transformation (34) has an 
inverse t(x), since dx/dt > 0 V t > t^. We also define 

Xq = x(t=0). (35) 

We observe that for t^ ^ 0 we have x^ ^ 0. Recalling the constant-coefficient results 
we will call the quantity /a(t)b(t) the '^intensity of combat** ( see also Taylor and 
Parry ^ ; since the larger it is, the more quickly the battle is moving towards ter- 

mination. The average intensity of combat is given by /a(t)b‘<t) = (1/t) | /a(s)b(s)ds. 
Then we have ^ 

X - Xq = {(1/t) I /aTsVbTs) ds}t = /a(t)lT(t) t. (36) 

0 

The substitution (34) transforms (2) into 

d^x/dx^ + (l/2){d 5,n R(t)/dx}dx/dx - X = 0, (37) 

with initial conditions 

1/2 

x(x=Xq) = Xq, and (R (t)dx/dx>^^^ = -y^. 

Theorem 1 tells us that unless (37) is a constant-coefficient equation, it is impossible 
to transform the X force-level equation (2) into a constant-coefficient equation by 
a transformation of the independent variable alone.. Also, equation (37) is highly 
significant because it clearly shows us that the course of combat depends on just two 
weapon-system parameters: (1) R(t) = b(t)/a(t), the relative fire effectiveness 

(X to Y) of the two combatants, and (2) I(t) = /a(t)b(t) , the intensity of combat 
(through equation (34), which relates I(t) to x) . Both these parameters may vary 
over time. In particular, from (37) we see that the nature of temporal variations in 
relative fire effectiveness will have a significant effect upon the course of combat. 

For the power attrition-rate coefficients with no offset (7) , the transformed 
X force-level equation (37) becomes 
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d^x/dT^ + { (2q-l)/T}dx/dT - X = 0, 



(38) 



with initial conditions 

x(t=Tq) = Xq, and ^dx/dx}^_^ = -y^ (/kak^/ (y+v+2) ) , 

where q = (v+1) / (p+v+2) and 

T = T(t) = ( 2 / ( m+v+ 2 ) ) ( t+Kg ) ^ ^ ^ . (39) 

Hence, Tq = ( 2 A^k^ / ( y +v+2 ) ) Kg ^ V*+'^+2 ) / 2 ^ observe that Vp,v > -1 we have 

0 < q < 1. Furthermore, q > 1/2 ^ dR/dt > 0, i.e. R(t) is a strictly increasing 
function of time. 

5. Lanchester-Clifford-Schlafli Functions 

Consider the function defined by the power series 

00 

F^(x) = r(a) I (x/2)^^{k!r(k-Hx)}. (AO) 

k=0 

For a ^ 0,-1, -2,... the radius of convergence for infinite by the ratio 

test for convergence of power series ( see , for example, KNOPP^^^^). Hence, 
is an entire function of the complex variable z = x+iy with an essential singularity 
at the point at infinity. Now consider the function defined by the infinite 

series 

CO 

H (x) = r(a) I (x/2)^^^‘’’“V{k!r(k+«+l)}. (41) 

k=0 

Observing that 

H^(x) = (1/a) (x/2)^“f^^j^(x) , (42) 

we see that for a > 0 the infinite series (41) is uniformly convergent on compact 
subsets of the complex plane. From (42) we can readily deduce the recursive relation 



F^(x) = F^^^(x) + {(x/2)2/[a(a+l)]}F^^2<^>- 

We will call the functions F (x) and H (x) Lanchester-Clifford-Schlafli functions 

a a 

( see Note 11). Other properties are readily deduced and are given in Table IV. 
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Table IV. Properties of the LCS Functions 

1. dF^/dx = (x/2)^ ^“h^(x) 

2. • dH^/dx = (x/2)^“"^F^(x) 

3. F^(x)Fj^_^(x) - H^(x)H'^_^(x) = 1 Vx 
where ot is not an integer (including zero) 

4. F (x=0) = 1 

5. H (x=0) = 0 for a > 0 

a ^ 

6. dF /dx(x=0) = 0 

a 

7. {(x/2)^“^“dH^/dx}^^Q = 1 

8. ~ cosh X 

9. ~ sinh x 
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The function satisfies the second order ordinary differential 

d^F /dx2+ {(2a-l)/x}dF /dx - F = 0, 
a a a 

with initial conditions 

F^(x=0) =1, and dF^/dx(x=0) = 0, 

while H (x) satisfies 
a 

d^H /dx^ - {(2a-l)/x}dH /dx - H = 0, 
a a a 

with initial conditions (for a > 0) 

H (x=0) = 0, and {(x/2)^"^“dH /dx} „ = 1 

a a x=0 

Thus, {F ,H- } is a fundamental system of solutions to 

’ a 1-a 

d2F/dx2 + {(2a-l)/x}dF/dx- F = 0, 

with Wronskian = (x/2) . Let us observe that (see Table III) 

Cj^(t) = F^(x(t)), S^(t) = { /k^k^/ (p+v+2) (T (t) ) , 

C^Ct) = Fp(x(t)), SyCt) = { ^kgk^/ (y+v+2) (T (t) ) , 

where p = 1-q. If we define 



T (x) = H- (x)/F (x) , 
a' ' 1-a' ' a ’ 



then 



Djj(t) = Tjj(t) = S^(t)/C^(t) = {/k^k^/ (u+v+2)}^*^ ^Hp(x(t) ) /F^ (x (t) ) 

where T^Ct) denotes a hyperbolic-like GLF. Observing that lim x(t) = -H», 
^ t->+” 

that is a strictly increasing function of x Vxe[0,+«) and 

0 & T^(x) < r(l-a)/r(a) for 0 ^ x < -H», 



with 



lim T (x) = r(l-a)/r(a) , 

x-»-+<o “ 



equation 

(44) 



(45) 



(46) 

(47) 

(48) 

(49) 

(50) 
we see 

(51) 

(52) 
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since by the value of Q* determined by Taylor and Comstock 



for the power attrition 



[ 27 ] 



rate coefficients (7), denoted as Q*(ij,v,Kq= 0) , we have ( see (30) and (50)) 
lim Tj^(t) = 1 /Q*(vj,v,Kq= 0) = { /k^k^^/ (y+v+2) (p) /F (q) . 

Comparing (38) and (46) , we see that the solution to (38) is given by 



x(t) . x„{Fp(T„)F^(T(t))-H^(T„)Hp(T(t))) 

,2q-l 



^ {F (Tq)H (x(t))-H (Tq)F (T(t))}. 



The time to annihilate X is determined by x(t=tY) = 0 and thus 

A. 

T^(t( 4» = {«„Fp(T„)+yo^7i^(.^/(u+v+2))‘'-''Hp(Tj,)}/ 

where is given by (49) and from (51) 

0 S t^(t) < r(p)/r(q). 

For tQ = -Kg = 0, (54) simplifies to 



T (T(t^)) = (Xq/Yo) (y+v+2) )P’‘^ . 



From (54) and (55) we may deduce the following theorem: 



THEOREM 2: Consider combat between two homogeneous forces described by 

(1) with power attrition-rate coefficients (7). Assume that these 
equations hold for all time and that Y **wins” when ^ with 

y(t^) > 0. Then Y will win if and only if 

r (q) {xgFp (TQ)+yQ^t“/k^(,/^/ (p+v+2) (Tq) } < 

r (P) (To)+yo'^ka^^j^’'^/ ) ‘^■Pf^ (Tq) } . 

For t =0 (i.e. K„ = 0 and x_ = 0) , Y will win if and only if 

0 U 

XQF(q) < yQ/k^/k^( /kakj^/ (y+v+2) )^^ Pf(p). 



(53) 



(54) 



(55) 



(56) 
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6. Tabulations of LCS Functions 



Tabulations of the Lanchester-Clif ford-Schlafli functions are available in two 

of the authors’ reports, also available from the National Technical Information Service 

( see references 25 and 26) . These reports contain five-decimal-place tables of the 

hyperbolic-like LCs functions F (x) , H (x) , and T (x) for values of the argu- 

ment x = 0.00(0.01)2.00(0.1)10.0 and various values of the order a. The short 

table^^^^ contains tabulations for a = 1/2 ,1/3 , 2/3 ,1/4 ,3/4 ,1/5 ,2/5 , 3/5 , 4/5 ,3/7 , and 

f 261 

4/7 corresponding to y,v = 0,1, 2, 3; while the longer table^ ^ contains tabulations 

for a = 1/2,1/3,2/3,1/4,3/4,1/5,2/5,3/5,4/5,2/7,3/7,4/7,5/7,4/9,5/9,3/11,5/11,6/11, 

8/11,5/13,8/13,5/17,12/17,5/21, and 16/21 corresponding to y,v = 0,1/4 ,1/2 ,1,1 1/2, 

2,3. As we have seen above in Section 1 ( see (5) and Figure 2), such values of y and 

V allow one to analyze, for example, a wide variety of range capabilities for weapon 

systems in Bonder ’ s constant-speed attack model (5). These tables have been 

[24] 

calculated by the recursive methods given in Section 8 of Taylor and Brown 

A representative tabulation of the hyperbolic-like LCS functions F^(x) , 

H^_^(x) , and ^ ” 3/5, similar to those that appear in references 25 and 

26, is given in Tables V and VI. The values of the argument x are the same as those 
used for the tabulation of the hyperbolic functions by ABRAMOWITZ and STEGUN^^^. We 
observe from Table VI and (52) that the limiting value of x (here 

ot = 3/5) is quickly reached, with three-decimal-place agr.e^ment by x = 4.5. 

7. Numerical Examples 

In this section we examine a couple of numerical examples to show some of the 
insights that may be gained into the dynamics of combat between two homogeneous forces 
from our new results. As in references 21 and 24, we consider S. Bonder model 
(5) for the constant-speed attack against a static defensive position. We will focus 
on the new results of this paper [in particular, the prediction of battle outcome from 
initial conditions without explicitly computing the force-level trajectories (cf. 
questions (Ql) and (Q4) of Section 1)]. From the input data given in Table VII, we 



23 








®0>m®^ 


r-o®iTk® 


0^>®-4® 


®fNi9>9i,^ 


r4N#00# 


mmrsj#o 


— ®®— 


®® — FM® 


0>IAA-«A® 


r- 




lAl/%«A«S-4 


COlAOtfNO> 

^—'OOsr 




o><^ ►-»/%-*> 


0<0<M»*— 4 


lA®— mi/> 




m— <^®rsJ 


® ^o>®® 


rvi 






^m^>-4iA 


(^m^--4u^ 






A4kA® — #' 


>-OAJlA® 


omiA®o 


m 




00-4-4^4 


># 


lAlA® ® 


<o^-^•®® 


®0>0^CKO 


00-4-4 — 


<Njr\»f\jmm 


m*^i4»f?km 


lAlAlAlA® 


® 




rgr^rsjrsjiM 


rsjr>jf>j(\ir4 


<MrMr4rsjf>i 


mmmmm 


mmmmm 




m 




-4<^(*tOr^ 




m®iA-4^ 




■#rgO^-n 




C^^-®F-0 


lAAIfSiiAO 


9iAI<^(^lA 


lA 






rNiOO-4it 








®'#>r 


lAO>lAm<# 


4)— 


fMor-®oj 


® 




>r^>ru%K\ 


<«^<•®(^o 


— mK^^-(^ 


-4'#<«9‘<M 


•A<^fM'00 


^9>m®m 


®m9>kA— 


^•i#0— ® 


<^J<^^.^A® 


AJ 




3S23S 


0(M>#>0» 


— 

ooooo 


fsj^^®— 


mi/>®om 


t/>h-0«\(lA 


r>-0<N<»A® 


Om®®— 


®®Q»AJIA 


® 




9*Ch»0^0* 


— — IM 






®IAIAIAIA 


®®®®^• 


A-#-A-®® 


® 


• • • • • 


• • • • • 


• • • • • 


• • • • • 

rMf^«wrg 


• • • • • 


• • • • • 

<M4M<Mcgrvj 


• c • • • 

AjrsjAjrgfNj 


• • • • • 

rg«vjfvjrgrg 


• • • • • 


• • • • • 


rsj 










h.h.^®« 


-4o>0‘>r<'j 


(#— 


Ch®— O'rsJ 


o<#mh-^* 


m#^CMiAiA 


<MiA<# — <# 


lAm®— 


O®— ®IA 




N09^»0 

oo^d^o 








®— ®iAm 

omi«>®— 


J2251 


mmiA®m 


Q®®om 

•#®m(^>t 


«iAm^® 
<^iA— r^im 


ommrsjm 


® 


00-4 -4 rsj 




OOJ'T ■£«> 


Aj<oo4ro* 


o®moA- 








iA«ri/\iAK\ 


— 


rw®0*OfNi 


P2222 


o— m«#(A 

•«®®« 


^-®o-4m 

®®<^o*<^ 


®®®»-4 

^0'<^0*O 


m®®®» 

ooooo 






















AiCgfMAiOg 


• 

AJ 






0-4 rsim^t 




0-4<Nim'# 


tf\4i^-®0* 


0^tS4tfS^ 




o— Aim® 


lA®r-®|/> 


o 


ooooo 


ooooo 


-4-4—4 — 


— — — 










®®®®® 


®®®®® 


lA 




® ®f*^®-4 


— 0®(«>P- 


<^<^®®<^J 


®rsjiAh>® 


(>0>®®iA 


mo®iom 


0®®IA<# 


mm® ®o* 


rNjvom— o 


O 






®m®o>o 


0^h-i»0 


<#00-4-4 


o®m— ® 


om<#iAiA 

9><omo^- 


<#— «<#0^ 


m®®(^(^ 


(^^-l/^AJ® 


m 






OO0^®«P 




0«^-lAm 


— ®®<#— 




f^m^iA— 


h-mo^mo 


® 




®®®«® 


®d‘0^o— 

®®«»o> 


rsim^iA® 


>^►••0*0 

<^0>9>9'0 


— <Nim^ 
OOOOO 


ooooo 


®o^^o— 

000—4 


— (Mfsjm® 


®IAIA®A- 


A. 


OOOOO 


ooooo 


ooooo 


OOOOO 


0000-4 


















00>A®(*t 


IA'«>®9>^ 


<M'#mO'» 




<#0>rNjtA® 


-4>r®rNi® 




rsj®<oF-o 


r». 




®!«>®>rO 


^•^r-4(^® 


^-^-^-®o 


m<oomo 


®m— 09> 

®irkra»iA 


o— <#r- — 


^-mo(^® 


o*Omf-i*g 


<^®lAlAr- 


O 




®^0^lA-4 


<«(SJ®m(^ 


K\— fomo 




mo»»->##>j 


C^^^A<^JO 


®^-®m<M 


OO^®^-® 






mi/\®®0 


— 


S— — — — 


h-O^OfNJ<# 


l/\^-9>0<^J 


#>0^0‘-4 


rsj'#®®o 


— miAh-o^ 


— <M® ®® 


o 




O'O'^O'O 


OOOOO 


— — rgrsj<Nj 


o^rsjrsjmm 


mm mm <# 




lAlAlAlAlA 


®®®®® 


A- 


OOOOO 

<SJO^r»>0‘ 


0000-4 


— 4 — — 


— -4 — — 


— — 


— M— — 


— — — — 


— — — 




— — *— * 


* 


IAOI/>0‘X^ 


►•o^rr-o 


m -oomf*- 


— mviAo 


® mo^® 


®(^ — #® 


mo®®(> 


iMf-m— 


m® ®m® 


® 


<M®O'0<V 


00>®®0 


<M<OOi/xrM 


o>r^p-r»® 


-4^®^ — 


®»^f>-r«-a> 


rsj<0<Ni®iA 


®®<#<MO 


® 0 ‘®®m 

*h-®iA® 


m®r^-4iO 


AJ 


'OOiAO^'T 








04®<#—® 


<#-4®m<M 


Or-iA<MO 


mrg— o 


o 








®(DO^O^O 


— — rgmm 


^^lAtrvo^ 


®® 9>0 — 


— fNjm#iA 


iA®r»®o^ 


0-4 Aim® 


lA 


• 1 lull 1 ^ 


>--4(14-4-4 


— — — 


— — rsj 


<Mr\j<Mrvr>j 


rvirs«M<MfM 


f>jrMrMmm 


mmmmm 


m i^m m#^ 


®®® ®-J 


® 






















— 


o-*FNjrn>r 




0-4fSJ|4>"#- 


l/^’0^-®(^ 


O— r>jm^ 




o— Ajm^t 


iA®r-® 


o— Aim® 


iA®r>-®9‘ 


o 


iAiAlAi/M/\ 


lAlAlAlAlA 


® 'O ® 


'O'O^'O'O 






® ®<u® ® 


®®®® ® 






o 


OOOOO 


OOOOO 


ooooo 


ooooo 


ooooo 


OOOOO 


ooooo 


ooooo 


OOOOO 


ooooo 






mo® ®*^ 


mo-4(^® 


^ 

o«rtr,<#o 


Om o— ^>- 


tA'OO^o® 


rv«rNjiA#® 


^-® 


®00A-0 


o®oo« 


AJ 


Oh~® rsj 


■«0(M® or- 


oa*-#*Ai# 


«rvj 


^0<VJ— V 


<o— it®® 


»nma>m® 




A*0*A®m 


AJ 


vOfvJ'OO' 


0-400® 






rsje-lvm-iM 


>0-4i/\0^fsj 


®om®9* 




®®®CrM 


mi/>®A»® 


0^ 


m®«o 

oooo-^ 




fvjit 

rsirg<MfNj(M 


— rsj4’ 0^- 
mm 


(^Of^JmtA 

m<#<#<##^ 


<0®^0<M 


miA®r^« 

lAiAiAiAiA 


O— 


®A-90— 


Aim®iA® 
r>r>.— h~r. 


Jk. 


4®®®® 


®®®#-A- 


f- 


ooooo 


OOOOO 


OOOOO 


ooooo 


boobo 


OOOOO 


OOOOO 


bboob 


OOOOO 


ooooo 


o 


►-o**»>* 


®mm®® 


®m® ^--4 


m®-4— 1/\ 


omrao*^- 


^-® rvjOrg 


<^m# 


rsih--4m® 


® ®mo4Aj 


miA®mo 




00<»f«> 




> — 0>IA<^ 


00*f-m^- 


0<Mmr4-4 




«m®om 


®®— miA 


iw(^— miA 


A>®>— ®r- 


O' 


O'MOO' 


0-4— oo> 


® 


^0<M9>'0<V 


^i/\— f»m 


®'tOl/\-4 


^r\ff-m® 


m®<#<^® 


9«®0«A0 


lAO® — ® 




m>o®0 






— m^<0® 
mmrnmm 


0>— m#® 


r»>o>— (M^ 


lAP-®0 — 


m®®^-o^ 


OAJ®iAA- 


®o-4m® 


® 


0000-4 


>4>— — ^ 


<NrsjfsjfNj<M 


**'>#'#<#'# 




lAlAlA®® 


®®®®® 




A>® s® ® 


« 


• • • • • 


• • • • • 


• • • • • 


• • • • • 


• • • • • 


• • • • • 


• • • • • 


• • • • • 


• • • • • 


• • • • • 


• 


ooooo 


OOOOO 


OOOOO 


OOOOO 


ooooo 


OOOOO 


ooooo 


ooooo 


ooooo 


ooooo 


o 




'ro■^■^-® 


r-tf>-4l/>® 


9>9>f^m® 


Aj.#«#m-4 




® <#— ®o 


miAiAtAm 


OA-<M®9> 


Aim® ®m 


AJ 


00-4 r0 4> 


OiAO>Of<^ 


-4000— 


m®OiAO 

O^OtsimiA 


Ki.#fSi-4 — 


— miA®<M 


r*-mQr-® 

h-omiA® 


lAlA®® — 


iA<MA— ® 




AJ 


OOOOO 


— — 




<090fM<# 


®®Orsji/\ 


— ^r^O'# 


^-0®®-4 


ia<a mr-— 


® 


ooooo 


OOOOO 


ooooo 

OOOOO 


Q— 4 — 


— — rsirvjrsj 


(Njrsimmm 


ooSoS 


lAiAlA®® 


«r-rofo0 


SSSS2 


o 


OOOOO 


OOOOO 


ooooo 


OOOOO 


OOOOO 


OOOOO 


OOOOO 


— 

• 


— <Mrci^ 


IA®^®0> 






o-4rgm<r 


iA®fk.®9> 


o-4ram>» 


IA®^>®9> 


0-4»M|4>it 


IA®A.®® 


o 


OOOOO 


OOOOO 




m $m $i t 


fsjrgrg«M^ 


o<<Njr4rM<sj 


mmmm m 


mmmmm 


®®®®® 


®®®®® 


lA 


• • • • • 


■ • • • • 


• • • • • 


• • • • • 


• • • • • 


• • • • • 


• • • * • 


• • • • • 




* • • • • 




ooooo 


OOOOO 


OOOOO 


OOOOO 


OOOOO 


OOOOO 


OOOOO 


OOOOO 


OOOOO 


OOOOO 


o 



24 



Table V. Lanchester-Clif ford-Schlaf 11 Functions F (x) , H, ^(x), and 
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Table VI, Lanchester-Clif ford-SchlSfli Functions F (x) , H. (x) , and 



Table VII. Input Data for Numerical Examples 



p = 1, V = 2 

Oq = 0.06 X casualties /minute /Y unit 



= 0.6 Y casualties/minute/X unit 



R = R_ = 2000 meters 



V = 5 miles/hour 
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compute the parameter values shown in Table VIII. We observe from Tables VI and VIII 
the predicted agreement between r(l-a)/r(a) and the limiting value of T^(x) as 
X -> +<» ( see (52)) for a = q = 3/5. We now consider two cases: (I) = 2000 

meters, and (II) Rq = 1250 meters. 

[211 

When Rq = 2000 meters ( see Figure 3 of Taylor^ O, we have = 0 and 

= 0. The maximum time that the battle can last is t =14.91 minutes, since at 
0 max 

this time the attackers reach their final objective (i.e. the defensive position). 

We now consider the qualitative behavior of the M= 1, v = 2 force-level trajectory 

[ 21 ] 

shown in Figure 3 of Taylor . Theorem 2 tells us that X can be annihilated 
^ ^ 0-^20. By (56) the annihilation time of the X force is given by 

Tq(x(t^)) = 3.544 ^'or Xq = 10, yQ = 30, we have T^(t^) = 1.18122 so that 

from Table V (using linear interpolation) we obtain T = 1.009. Hence, (39) yields 

A 

a a 

t^ = 14.24 minutes and r^ = 89.8 meters. Further results are given in Table IX. 

[241 

When Rq = 1250 meters ( see Figure 3 of Taylor and Brown ), we have 

= 5.5923 minutes, T = 0.0975, and t = 9.32 minutes. Jin this case (again, 
o U max 

for y = 1, V = 2), X can be annihilated ^ 0*382 with [from (5A) ] the 

annihilation time of the X force given by T (t^) = (3.565 u +0.223)/(0.156 u +1.004), 

q A u u 

where Uq = ^q/Vq^ Some further numerical results are given in Table X. Again, these 
parametric results should be contrasted with the single y = 1, v = 2 force-level 
trajectory shown in Figure 3 of reference 24. 

8. Discussion 

In Section 7 above we have seen how our new definition of power Lanchester 
functions (guided by the general requirements for GLF given in Table I) allows one to 
conveniently obtain much valuable information about the model (1) with attrition-rate 
coefficients (7) without explicitly computing the entire force-level trajectories ( see 
Note 12). Previously we were limited to only computing force-level trajectories. 

Now we can tell who is going to be annihilated and when without explicitly computing 
the trajectories ( see Note 13). Not only did we answer questions about qualitative 
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Table VIII. Parameter Values for Numerical Examples 



k = 4.0233 x 10 ^ X casualties/ (minute) ^/Y unit 
a 

= 2.6979 x 10 ^ Y casualties/ (minute) '^/X unit 
p = 2/5, q = 3/5 
r(p)/r(q) = 1.48951 
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Table IX. Annihilation of the X Force as a Function 
of the Initial Force Ratio for Rq = 2000 meters 





t^ (minutes) 


r^(meters) 


0.333 


14. 2A 


89.8 


0.250 


11.61 


443.2 


0.200 


10.19 


633.2 
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Table X. Annihilation of the X Force as a Function 



of the Initial Force Ratio for Rq = 1250 meters 





tjj (minute's) 


r^(raeters) 


0.333 


10.63 


t 


0.250 


7.56 


235.9 


0.200 


6.17 


422.8 



t =9.32 minutes 
max 



and = x(r=0) = 1.35. 
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model behavior (e,g. force annihilation) for specific values of, for example, initial 
force levels but also for a range of values of the initial force ratio (i,e, parametric 
analysis of model behavior). The results of this paper may be used for other para- 
metric analyses ( see Note 14), e,g, parametric dependence of battle outcome on attri- 

[241 

tion-rate coefficients. Thus, our extension of past results allows one to develop 

important insights into the dynamics of combat between two homogeneous forces with 

r 2 s 1 

temporal variations in fire effectivenesses. With the availability * of tabula- 
tions of the LCS functions, one can now analyze such combat modelled by the power 
attrition-rate coefficients (7) with somewhat' the same facility as he can for the 
constant-coefficient case ( see Note 15) and thus aid in parametric analyses. 

In his classic 1914 paper Lanchester assumed that the combatants* fire 
effectivenesses (as expressed by Lanchester attrition-rate coefficients) were constant 
over time and deduced his famous square law 

e{x2-x2(t)} = a{y2-y2(t)}, (57) 

where a and 6 denote constant attrition-rate coefficients. It follows from (57) 
that (provided there is no **time limit** for the battle) 

X will be annihilated ^ 

Thus, we see that equality of Lanchester-type fighting strengths depends on two param- 
eters: (I) initial force ratio, and (II) relative effectiveness. When the timing 

of military actions is considered, we add a third parameter, the intensity of combat = 
to this list of significant combat parameters. No such simple relationship like 
the square law (57), which yields (58), holds in general foi* variable attrition-rate 
coefficients when a(t)/b(t) ^ constant. However, by transforming the independent 
variable t to normalize the battle’s time scale by the intensity of combat, we found 
( see equation (37)) that the course of such variable-coefficient combat depends on 
only two weapon-system parameters: (I) relative fire effectiveness, R(t) = b(t)/a(t), 

and (II) intensity of combat, I(t) = /a(t)b(t) , Moreover, we extended (58) to combat 
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modelled with the power attrition-rate coefficients with **no offset’* (7) ( see Theorem 2). 
This is the first time that such a generalization of the square law has been obtained 
for the variable-coefficient Lanchester-type model (1) with a(t)/b(t) ^ constant. We 
observe that for > 0 this ’’exact” outcome-prediction relation (i.e. necessary 
and sufficient condition for force annihilation) involves higher transcendental functions 
(here, the LCS functions) and is complementary to the sufficient condition (involving 

r 23 1 

only elementary functions) given by Taylor and Parry for > 0. 

F5 71 ri21 F2 91 

Work by BONDER * , Clark^ , and others*’ * ^ on the prediction of Lanchester 

F241 

attrition-rate coefficients ( see Taylor and Btown for further discussion and 
references) has generated interest in variable-coefficient Lanchester-type models. 
Interest in the power attrition-rate coefficients with ”no offset” (7) is provided by 
S. Bonder ’ s model (5) and his examination of predicted attrition-rate for various 
weapon systems ( see pp. 196-200 of reference 9). However useful our results may be 
in their own right, they have far greater import: (I) they are a model for the treat- 
ment of other Lanchester functions and their tabulations, and (II) they may be used in 

F27 1 

the numerical determination of the parity-condition parameter Q* for related 

attrition-rate coefficients (for example, (4) with > 0) . In the future we will 
show how our tabulations of the LCS functions play a key role in the numerical deter- 
mination of the parity-condition parameter Q* for the general power attrition-rate 

coefficients (4) with positive ’’offset” (i.e. Kq > 0) . 

F241 

We have extended our mathematical theory of variable-coefficient Lanchester- 
type equations of ’’modern warfare” for combat between two homogeneous forces in order 
to be able to more thoroughly analyze such models. The classic ordinary differential 
equation theories ( see , for example, HILLE^^^^) were inadequate to supply all the 
answers sought about such combat models (cf. questions (Q1)-(Q4) in Section 1 above). 

The mathematical theory of the model (1) with coefficients (7) is now nearly as complete 
as that of the constant-coefficient model. Such results as we have given here are very 
useful for understanding the dynamics of combat, i.e. how the trading of casualties 
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will be projected over time. H. K. WEISS has emphasized that such a simplified 
model of a combat situation is particularly valuable when it leads to a clearer 
understanding of significant relationships that would tend to be obscured in a more 
complex model. As is always the case, however, the insights geined into combat 
dynamics are no more valid than the models themselves. 

9 . Summary 

In this paper we have introduced new mathematical functions (Lanchester-Cliff ord- 
Schlafli, or LCS, functions) that allow important information (in particular, force- 
annihilation prediction) to be obtained without explicitly computing force-level 
trajectories for the variable-coefficient Lanchester-type model (1) with power attri- 
tion-rate coefficients with "no offset" (7). Our development was based on new theoretical 
considerations: we gave a new general discussion of representing the solutions to the 

X and Y force-level equations in terms of general Lanchester functions (GLF) and 
gave the general properties that these GLF should possess; we shox^ed that there are 
essentially only two kinds of GLF that satisfy these requirements (exponential-like 
GLF and hyperbolic-like GLF) and that the hyperbolic-like functions are to be preferred. 
Moreover, the exponential-like GLF are an essential theoretical construct, since they 
play a key role in determining force-annihilation-prediction conditions (i.e. showing 
that the reciprocal of the parity condition parameter is equal to the limiting value 
of the quotient of two hyperbolic-like general Lanchester X-functions) . We stressed 
that such building blocks should be chosen to yield as much information as possible 
about the model (and as conveniently as possible). We saw that the analysis of, for 
example, the X force-level equation was facilitated by transforming the battle’s 
time scale and that the only two weapon-system parameters affecting the course of 
combat are the relative fire effectiveness and the intensity of combat. These results 
extended and unified our mathematical theory of variable-coefficient Lanchester-type 
equations of "modern warfare" ( see reference 24) . 
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We then applied our general mathematical theory to the special case of combat 
modelled by power attrition-rate coefficients with ”no offset.” Our new definition 



of Lanchester-Clif ford-Schlaf li (LCS) functions was required for these power attrition 
rate coefficients in order to answer questions about battle outcome without explicitly 
computing force-level trajectories (i.e. to predict battle outcome/force annihilation) 
The mathematical theory of this variable-coefficient Lanchester-type models of "modern 
warfare” (modelling, for example, weapon systems with the same effective range) is 
now nearly as complete as that of the constant-coefficient model. With tabulations 
of the new LCS functions now available, one can study this variable-coefficient model 
almost as easily and thoroughly as Lanchester’s classic constant-coefficient model. 



1. Following terminology introduced in reference 24, we will refer to Lanchester 
functions corresponding to the power attrition-rate coefficients (4) with Kq > 0 as 
offset power Lanchester functions ( see Section 1). The power Lanchester (i.e. LCS) 



"no offset” (7). 

2. The equations (1) are only valid for x,y > 0. The first, for example, becomes 
dx/dt = 0 for X = 0. 

3. Further information on sets of circumstances that have been hypothesized to yield 



Notes 




the combat equations (1) (with constant coefficients) may be found in BRACKNEY 



[ 11 ] 



and Weiss 



4. It is impossible for all a(t),b(t) > 0 satisfying Condition (A) to have 



x^,y^ > 0 for i = 1,2 and all t > t^ such that 




and “ 




b(t)x^ and x^,X 2 are linearly independent. If it were possible, then 
|S(t)| = without loss of generality we may take Cg > 0. 
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Introducing and A = . we would have |S(t)l = y-^2^ “ *"S ^ ® 

A = /kj^/k^ b(t) (uj^+U 2 )A > /k^/k^ b(t)A^ > 0 for t > so that A(t) > 0 is strictly 



Increasing for t ^ t^. It would follow that l/ACt^) - 1/A(t) ^ 



which is impossible for b(t) such that lim 

t->+oo 



b(s)ds = +®. 






b(s)ds. 



5. To keep X 2 > y^^, and as "simple" as possible we specify that their initial 

values at t^ to be either 0 or 1. In order that |s(t)| = ^ 1^2 ” ^2^1 

must therefore have either Xj^(t=tQ> = 72 ~ ® X 2 (t=tQ) = yj^(t=tQ) =1, or 

Xji^Ct-tg) = y 2 (t=tQ> = 1 and X 2 (t=tQ> = yj^Crpt^) = 0. We consider the first possi- 
bility with similar arguments holding for the second. In this first case examination 
of the differential equations with initial conditions shows us that x. = S , x = C , 
~ ^2 ” functions coincide with the hyperbolic-like GLF of 

Table III. Thus, we need not consider L = I, since the same results may be obtained 
by using another one of Lemma I’s feasible forms for L. 

6. We conjecture that some condition like lim T(t) = -h» is sufficient to guarantee 

t-H-“ 

that Q* is unique. 



7. This intuitively obvious result may be proved by observing the identity 

rt 



{b(o)x2(a) + a(o)y2(o)}do = x(s)y(s) - x(t)y (t) . 



a(t)dt and 



^0 



b(t)dt is 
[15], 



A less obvious fact is that unless at least one of 

[15] 

unbounded, then neither x(t) nor y(t) need ever be annihilated (see Hille ). 

As an example of this situation, consider the batt.le with attrition-rate coefficients (7), 
Kg > 0, and p = v < -1. Then x(t) = Xq cosh 6 (t) - yQ^k^T^sinh 6 (t) , where 6(t) = 
{-l/Cv+DHl/K^"^^'*'^^ - l/(t+K_)"^'^'*'^^ . Let e = {-l/(v+l)Hl/K„“^'''‘’^h > 0 and 
finite. It follows that v < -1 and Kg > 0 can be chosen so that even though 



Xq < y^v^^/k^ we have lim x(t) = x^ cosh 0^ - y^/k^Tl^sinh 0^ > 0. 



35 



8. In general, the value of Q* will not be known exactly. Unfortunately, errors 



in the initial conditions for the exponential-like GLF become exponentially magnified 
over time. The situation is even worse for Ti^(t;Q*) = E 2 ^(t;Q^) = E^(t ;Q*)/E^(t ;Q*) , 
which is used to determine the time that the X force will be annihilated. 



[221 [231 

9. The term quasi-autonomous was coined by Taylor ( see also TAYLOR ) to denote 
a system of differential equations that may be transformed to an autonomous system ( see 

ri8i 

for example, p. 163 of PETROVSKI^ by a change of the time scale. Special cases 

[9] 

of such Lanchester-type equations have been considered by, for example, Farrell and 
TAYLOR^^^^. More general (possibly nonlinear) quasi-autonomous Lanchester-type equa- 
tions have been studied by Taylor ( see also Note 4 of Taylor and Brown^^"^^). 



10. If we were to specify the initial conditions of the GLF at t = 0 instead of 

t = t^, then (32) would reduce to x(t) = XQC^(t) - * However, when the 

initial conditions for the hyperbolic-like GLF are not given at t^, a separate tabu- 
lation of, for example, must be used for each different value of t^ 

(i.e. = C^(t5tQ)). 

11. Although the solution to the X force-level equation (2) with the power attrition- 
rate coefficients (7) may be expressed in terms of known higher transcendental functions 
(see Taylor , Taylor and Brown^^^^, and Taylor and Comstock^^^^ ) , we have chosen 

to introduce the LCS functions, since tabulations of these other functions are not 
available for the full range of parameter values of interest in Lanchester combat 
theory. For example, we can construct such solutions with modified Bessel functions 
of the first kind of fractional order, but tabulations of these ( see , for example, 
Abramowitz and Stegun^^^) only exist for a restrictive set of values of the order p 
(i.e. p = ±1/4, ±1/3, ±1/2, ±2/3, ±3/4), where p = (y+1)/ (y+v+2) . Furthermore, 
tabulations of functions corresponding to the quotient of, for example, two GLXF do 
not apparently exist at all. Consequently, we have introduced our new LCS functions, 
which provide much of the information desired about such battles. The naming of our 
LCS functions follows from the facts that a function similar to F (x) was introduced 
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ri9i 

by LUDWIG SCHLAFLI (1814-1895) in 1867 > while a related one appears in a posthumous 
fragment of the great English geometer William Kingdon Clifford (1845-1879) ( see pp, 343- 
348 of CLIFFORD . 

[141 

12. In his well-known survey paper on the Lanchester theory of combat, Dolansky 
suggested the development of outcome predicting relations without solving in detail 
and/or computing force-level trajectories as one of several problems for future research. 
Our Theorem 2 is a step towards this problem's resolution ( see also references 22, 27, 
and 28). 

13. Bonder and Honig^^^^ point out, however, that force annihilation may not be the 

appropriate criterion for evaluating many military operations, especially when force 

[9] 

annihilation does not occur. See pp. 192-242 of Bonder and Farrell for a detailed 
Lanchester-type analysis of an attack situation for which other ”end of battle condi- 
tions** play the major role in the evaluation process. Nevertheless, it is of interest 
to know when and why force annihilation will occur. 

[81 

14. S. BONDER^ ^ has suggested that an increased emphasis be placed on parametric 
analyses in systems analysis studies ( see pp. 21-22 of reference 8). 

15. One significant exception is that the outcome of fixed-force-level-breakpoint 

battles (for example, Y "wins** when x^ = x(t^) = x^p but y^ > 7gp> where t^, 

Xf, y^ denote final values and x^p denotes X*s breakpoint) with x^p y^p > 0 and 

a(t)/b(t) ^ constant cannot apparently be analyzed in the manner described in this 

[27] 

paper (see Taylor and Comstock ). 
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